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ABSTRACT 


In  measurement  of  dissolved  oxygen  (DO)  stratifi- 
cations^ (The  lowest  point  is  chosen  usually  at  a 
depth  of  0.5  m above  the  sediment  surface.  DO 
concentrai.ions  in  the  overlying  water  just  above 
the  sediment  are  rarely  measured.  Therefore,  the 
DO  behavior  in  the  sediment-water  interface  is 
not  well  known.  In  order  to  investigate  DO 
behavior,  column  tests  on  DC  were  performed  by 
use  of  special  test  apparatus.  This  paper  deals 
with  these  test  results  and  discusses  the 
boundary  conditions  of  sediment-water  interface 
on  DO  behavior,  sf. 


OUTLINE  OF  TESTS 


The  test  column  consists  of  two  parts:  the  upper  and  lower  sections. 
The  latter  is  the  column  sampler  itself,  which  has  been  used  for  sampling 
undisturbed  sediment.  Both  tubes  are  connected  by  a ring  coupling  in  the 
laboratory  (Figures  1 and  2). 


The  connected  test  column  is  then  provided  with  water  jackets  at  its 
circumference  and  erected  in  an  isothermal  tank.  In  order  to  emerge  from  a 
thermal  stratification  in  summer  lake  water,  the  overlying  water  column  is 
heated  by  three  isothermal  tanks  with  different  temperatures.  Since  we  had 
to  measure  DO  concentrations  of  very  thin  layers,  we  have  paid  close  attention 
to  water  sampling  and  analysis. 
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Figure  1.  Test  apparatus 

The  sampling  tube  (Figure  3)  has  a diameter  of  6 mm,  a horizontal  length 
of  170  mm,  and  8 holes  of  0.5  mm<J>.  In  the  1980  tests,  4 tubes  were  inserted 
at  heights  of  1,  10,  40,  and  70  cm  above  sediment,  and  in  the  1981  tests, 

6 tubes  at  heights  of  1,  4,  7,  10,  40,  and  70  cm,  respectively.  The  15-ml 
water  samples  were  taken  gently  so  that  the  Reynolds'  numbers  did  not  exceed 
100.  Samples  were  put  into  the  special  DO  meter,  which  was  designed  for  this 
test  (Figure  4).  This  meter  has  a small  stirrer  for  medical  use  and  is  capable 
of  measuring  very  small  quantities  of  water. 

The  sediment  samples  were  taken  from  Lake  Nakanoumi  and  one  of  its  in- 
flowing rivers. 

THERMAL  STRATIFICATION 

Quantities  of  heat  flowing  into  and  out  of  a segment  of  the  water  column 
may  be  expressed  as  shown  in  Figure  5. 
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Figure  2.  Photo  of  apparatus 


Electrode 


Figure  3.  Sampling  tube 

As  equilibrium  equation  for  the  heat  quantities  is  represented  as  follows: 


where 


3T  k 32T  _ 3p  9T 

3-  = — 3 - wT  37  - pw  — 

9t  up  gz2  OZ  dz 


(1) 


T = water  temperature  °C 
t * time  after  incubation  of  tests,  days 
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Figure  5.  Induction  of  thermal  stratification 

w = descending  flow  velocity  due  to  water  temperature  difference,  cm/sec 

3 

Cp  = specific  heat  of  water,  cal/°C/cm 
k = thermal  conductivity,  cal/cm/sec/°C 
p = water  density,  g/cm^ 
z = height  of  water  column 

The  solution  of  Equation  1 is  obtained  in  the  following  series  function: 


T_ 

T 

o 


(2) 
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where 


Tq  =>  water  temperature  at  z = 0 

p20  = water  density  at  20°C 

p ■ P2Q  (1  - aT)  is  assumed 

«p20cp 
6 *- — 

2 

D “ thermal  diffusivity,  cm  /sec 
Here  we  calculate 


2.55  x IQ  ^ x Q.9982  * 1 
1.42 


1.793  x 10 


*C  cm 


(3) 

(4) 


The  observed  water  temperatures  are  distributed  as  shown  in  Figure  6. 
These  are  represented  in  Figure  7 in  terms  of  the  ratio  of  T/Tq. 


Hater  Taap. 


Figure  6.  Water  temperature  distribution 
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Putting  the  test  results  T/Tq  - 1.092  at  z * 40  cm  (1981)  into  Equation  2, 
we  can  calculate  the  values  of  w: 

2 2 3 

1.092  ~ 1 + “2^4“  " 2A  d‘793  * 10“4)  x w3  x 21. 73  x AO3 


0.092  - 400  ~ - 1.0995  x IQ-5  w3 


Neglecting  the  second  term  in  the  right  side,  we  get 

w2  _ 0.092  m , in-4 
“2  ^00^  2*3  * 10 


| - 1.517  x lo”2  cm"1 

Similarly  for  T/Tq  - 1.1085  at  z - 30  cm  (1980) 

J - 2.195  x io'2  cm"1 
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—3  2 

Taking  D = 1.5  * 10  cm  /sec,  we  get 

w = 1.517  x lO-2  x 1.5  x io-3  = 2.276  x i(f5  cm/sec  (1981) 

w = 2.195  x io-2  x 1.5  x io-3  = 3.287  x 10~5  cm/sec  (1980) 

From  the  water  temperature  profiles  it  is  found  that  the  descending 
velocity  (w)  due  to  temperature  differences  at  greater  differences  is  about  one 
and  a half  times  as  large  as  at  lower  differences.  (This  means  that  the  value 
of  w is  proportional  to  temperature  difference.  The  larger  temperature  differ- 
ence, the  bigger  w.) 

VARIANCE  OF  DO  AGAINST  TIME 

As  indicated  earlier,  sediment  samples  were  taken  from  Lake  Nakanoumi  and 
one  of  its  inflowing  rivers.  .Tap  water  was  used  as  overlying  water. 

The  relationships  .between  DO  concentration  profiles  and  time  transfer  are 
shown  in  Figures  8 and  9. 


co 


Figure  8.  DO  stratification 
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Figure  9.  DO  stratification 

As  the  figures  show,  the  initial  DO  profile  is  a nearly  vertical  line, 
located  almost  entirely  in  the  portion  overlying  the  sediment  10  cm  above  the 
interface.  DO  concentrations  decrease  rapidly  and  approach  within  10  cm  above 
the  interface.  We  can  see  clearly  that  DO  concentrations  decrease  gradually 
according  to  time  and  reach  a final  distribution  shifted  to  the  left  in  these 
figures.  Their  reduction  is  strong  at  the  lower  layers  and  weak  in  the  upper 
layers  (Figure  10). 


DO  STRATIFICATION 

From  our  test  results,  we  have  learned  that  we  should  consider  DO  strati- 
fication in  two  parts.  The  first  part  is  the  distribution  due  to  thermal 
stratification  occurring  in  the  portion  of  the  water  column  above  the  water 
layer  immediately  overlying  the  sediment.  The  second  part  is  the  distribution 
due  to  sediment  oxygen  demand  and  occurring  in  the  layer  immediately  overlying 
the  sediment. 
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DO  DISTRIBUTION  DUE  TO  THERMAL  STRATIFICATION 


An  equilibrium  equation  of  disolved  oxygen  is  represented  as  follows: 


where 


(5) 


Cj  = DO  concentration,  mg/£ 

w = descending  velocity  due  to  thermal  stratification,  cm/sec 

K “ diffusion  coefficient  of  dissolved  oxygen  in  the  vertical  direction, 

2, 

cm  / sec 

z = length  of  water  column,  cm 
t ■ time,  sec 


We  assume  that  the  solution  of  Equation  5 is  a separate  function  of 
temperature  (T)  at  time  (t)  and  as  follows: 
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Then  putting 


We  get 


C]  = T(t)  S(z) 


w 3S  9^S 

T^.I-  3z2 


S(z) 


-n 


T'  = -n  T 


K£|.w||  + n2s 
3z2  32 


= 0 


From  Equation  7 


2 

T = Ce"n 

Equation  8 is  rewritten  as  follows: 


t 


2 2 
3S  W3S.nl 

, 2 “ K 3z  K 

o Z 


(»*-!»♦£)• 

2 W n2 
m -*■  + *-• 


0 


0 


0 


(6) 


(7) 

(8) 


(9) 

(10) 


Thereupon 

S(Z)  = Cje®!2  + C2em2Z 
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Finally 


Cj  = e 


-n2t 


/ mlS 

\Cle 


m->Z 


+ C2e 


(ID 


At  the  water  surface  (z  = H where  H = total  length  of  water  column),  the 
following  condition  is  composed  of 


"C  - Kf  - 0 


2 / m^H  m2H\  2 / m^H  m2H\ 

we-n  t (Cje  + C2e  j - Ke"n  1 ^C^e  + C2m2e  1 


/ \ mlH  / \ m2H 

(w  + mjK  )e  Cj  + I W + m2K  le 


C,  = 0 


w + n^K  H 


'2  w + m2K 


(12) 


Then  we  get 


-n2t  . 
c,  = e C 


[mjZ  w + m^K  fmj-m2  1 a + n»2z  | 

6 + w + m2K  6 J 


This  is  an  equation  on  DO  variance  in  an  unsteady  flow. 


(13) 


In  a steady  flow,  the  following  solution  is  obtained  from  the  equation: 


3c, 


32c, 


w dz  K dz2 


(14) 


The  curve  of  Equation  13  or  14  is  connected  to  the  curve  of  the 
deoxygenating  zone,  which  is  dealt  with  in  the  next  section.  Consequently, 
the  arbitrary  constant  Cj  is  determined  by  it.  Then, 
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z=h 


= I C2 1 


z=h 


Cl 


K 

w 


c0 


sin  a 


(15) 


C 


1 


Co 


sin  a 


(16) 


DISTRIBUTION  DUE  TO  THE  DEOXYGENATION  OF  SEDIMENT 

The  test  results  showed  that  the  water  temperature  in  the  water  just  above 
the  sediment  (z  <_  10  cm)  is  constant  (Figure  11).  This  means  that  there  is  no 
descending  velocity  due  to  thermal  stratification.  The  variance  of  DO  concen- 
trations at  this  portion  is  represented  as  follows: 


Figure  11.  DO  distribution  in  deoxygenating 
zone,  c^  = DO  distribution  due  to  thermal 
stratification,  c2  = DO  distribution  due  to 
sediment  02  demand 


9c2 

3t~ 


Xc2 


(17) 


where 

c2  = DO  concentration,  mg /l 

2 

K = diffusion  coefficient  of  dissolved  oxygen,  cm  /sec 
X = deoxygenating  rate,  sec 
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The  solution  of  Equation  18  is  obtained  as  follows: 


since  C = C0  at  t = 0,  z = 0,  we  get 


Co 

sin  a 


Equation  18  becomes 


C2 


Co 


sin  a 


(19) 


This  is  an  equation  of  DO  variance  at  the  deoxygenating  zone  in  an  un- 
steady flow.  In  a steady  flow,  it  becomes  as  follows: 


(20) 


X 


11.5  ~ 0.2 

11.5  x 3 x 24  x 3600 


3.73  x IQ-5  sec"1 


Putting  the  data  together: 
z = 70  cm;  Cj  = 8.8  rag/i  (t  = 0) 
z = 10  cm;  Cj  = 7.6  mg/l  (t  = 0) 
We  get 


(21) 


| «=  0.00244 
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Since  w is  already  known  (page  7),  then  we  obtain 

K = 1.345  x 10~2  cm2 /sec 

From  the  observed  data, 

t = 0,  C2  = 11.5  mg/i,  and  t = 3 days;  C2  = 0.2  mg/*,, 

we  get 

yj^=  0.0527  cm/sec 
aQ  = 0.025  rad 

Cq  = 1.22  mg/£ 

Figure  12  shows  the  DO  profile  in  the  deoxygenating  zone. 


L R 


DO 

Figure  12.  DO  profile  in  deoxygenating  zone 

BOUNDARY  CONDITIONS  ON  DO  AT  SEDIMENT  SURFACE 

As  the  test  results  show,  we  could  reemerge  DO  stratifications  in  the 
laboratory.  These  behaviors  were  analyzed  theoretically  in  the  preceding 
chapter. 
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Owing  to  the  DO  stratification,  the  diffusion  of  dissolved  oxygen  travels 
from  the  upper  layers  to  the  lower  ones,  that  is,  from  higher  concentrations 
to  lower  ones.  The  quantities  which  reach  from  the  top  to  the  sediment  sur- 
face are  -K  ^ per  unit  area  and  per  unit  time.  On  the  other  side,  some  quan- 
3z 

tities  of  oxygen  are  consumed  by  sedimeut  during  deoxygenation,  as  far  as  it  is 
polluted. 


These  two  quantities  should  be  in  equilibrium  (Figure  13).  Then  we  get 
the  equation  if  we  represent  sediment  oxygen  demand  (SOD)  as  0^. 


Figure  13.  Equilibrium  at  sediment  surface 


0 


2 
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(■ 


X - n 

cos  — £ — z + a 
sin  a 


■) 


0 


- n2) 


cot  a 


.-A 


(23) 


This  is  a fundamental  equation  of  SOD.  From  this,  we  can  say  that  SOD 
behaves  exponentially  against  time  (Figure  14). 
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Figure  14.  Sediment  oxygen  demand 
behavior  against  time 

The  total  consumption  of  oxygen  (Figure  15)  is  obtained  by  a planimeter. 
If  we  represent  it  mathematically,  it  becomes 


Figure  15.  Total  consumption  of  oxygen 
Q = AH  [/c(t  = o)  dz  - /c(t)  dzJ 


In  Table  1,  the  measured  and  calculated  values  of  SOD  in  one  sample  are 
indicated.  The  calculated  values  are  obtained  from  the  assumption  that 

n2  = 0.11. 


Table  1.  Measured  and  calculated  values  of  SOD 


SOD 

2 

o2/o2(o) 

_t 

mg/m  /d 

Measured 

Calculated 

0 

400 

1.0 

1.0 

2 

320 

0.8 

0.803 

4 

250 

0.625 

0.644 

6 

210 

0.525 

0.517 

8 

150 

0.375 

0.375 
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Figure  16.  Comparison  of  measured  and 
computed  values  of  SOD 

We  can  see  from  Figure  16  that  there  is  a good  agreement  between  observed 
and  calculated  values  of  SOD. 


CONCLUSIONS 

The  DO  stratifications  that  are  often  observed  in  lakes  and  other  water 
areas  are  well  known  to  pollution  investigators.  However,  there  are  very  few 
data  regarding  DO  behavior  on  the  sediment  surface.  From  our  test  results,  we 
can  say  that  the  sediment  surface  is  always  subject  to  an  anoxic  condition 
. owing  to  its  oxygen  demand. 

From  our  tests,  we  have  learned  that  any  DO  stratification  consists  of 
two  parts,  that  is,  one  due  to  thermal  stratification  and  one  due  to  the 
oxygen  demand  of  the  sediment.  It  was  also  found  that  the  latter  has  a small 
height  (z  £ 10  cm). 

In  winter,  no  thermal  stratification  occurs.  Consequently,  the  DO  profile 
becomes  a vertical  straight  line  in  most  of  the  water  column.  In  this  case,  it 
can  incorrectly  be  considered  to  extend  to  the  sediment  surface. 

The  line  bends  at  a large  angle  at  the  junction  with  the  deoxygenating 
zone  (Figure  17). 

DO  stratification  generated  in  summer  reveals  an  exponential  profile 


represented  by 


w 

Cj  — 

1 w 
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in  autumn 


Figure  17.  DO  profile 


This  curve  is  pulled  to  lower  concentrations  (to  the  left  in  the  figure) 
at  its  lower  end  by  the  curve  of  the  deoxygenating  zone. 

In  our  tests,  we  studied  SOD  behavior  from  both  observation  and  theoreti- 
cal considerations.  We  found  that  such  DO  behavior  on  the  boundary  layer  of 
the  sediment  has  a great  effect  upon  the  release  mechanism  of  nutrients.  How- 
ever, this  problem  was  not  within  the  scope  of  this  paper. 
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